High-energy Cosmic Rays and Neutrinos from Semi-relativistic Hypernovae 
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The origin of the ultrahigh-energy (UHE) cosmic rays (CRs) from the second knee (~ 6 x 10 17 eV) 
above in the CR spectrum is still unknown. Recently, there has been growing evidence that a 
peculiar type of supernovae, called hypernovae, are associated with sub-energetic gamma-ray bursts 
(GRBs), such as SN1998bw/GRB980425 and SN20031w/GRB031203. Such hypernovae appear to 
have high (up to mildly relativistic) velocity ejecta, which may be linked to the sub-energetic GRBs. 
Assuming a continuous distribution of the kinetic energy of the hypernova ejecta as a function of 
its velocity Et oc (Tf3)~ a with a ~ 2, we find that 1) the external shock wave produced by the 
high velocity ejecta of a hypernova can accelerate protons up to energies as high as 10 19 eV; 2) the 
cosmological hypernova rate is sufficient to account for the energy flux above the second knee; and 
3) the steeper spectrum of CRs at these energies can arise in these sources. In addition, hypernovae 
would also give rise to a faint diffuse UHE neutrino flux, due to pj interactions of the UHE CRs 
with hypernova optical-UV photons. 

PACS numbers: 98.70.Sa, 97.60.Bw 98.70.Rz, 
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There is a general consensus that galactic supernova 
remnants (SNRs) are responsible for the CRs at energies 
below the "knee" at ~ 3 x 10 15 eV most probably 
through the diffusive shock acceleration mechanism 0. 
Galactic SNRs expanding into their former stellar wind 
have been suggested to be responsible for CRs above the 
knee [![. Recent data from the KASCADE experiment 
suggest that heavy elements of nuclear charge Ze are 
accelerated by the galactic SNRs to the magnetic rigidity 
limit - 3 x 10 15 Z eV jj. Thus galactic SNRs may be 
able to produce CRs to at least ~ 10 17 eV (see e.g. [a, Q 
for recent reviews), an energy slightly below the "second 
knee" in the CR spectrum at ~ 6 x 10 17 eV (see e.g. Q). 
On the other hand, the highest energy CRs above a few 
times 10 19 eV are generally thought to be extra-galactic 
in origin, due to their isotropic distribution and a lack of 
galactic source candidates capable of producing them. At 
these energies, the possible sources include cosmological 
GRBs [§L a ctive galactic nuclei 0] and powerful radio 
galaxies [lfj |. 

The origin of the intermediate energy range CRs, 
10 i7 _ 1Q i9 eV) noweverj 

remains more elusive. Some 
authors have suggested a common galactic origin for all 
CRs between the first knee and 10 20 eV in young neu- 
tron stars or Magnetars , while others favor an extra- 
galactic origin of all CRs above 10 18 eV @, i, EH. Re- 
cent HiRes data shows that the transition from heavy 
nuclei to the proton composition may already occur at 
the second knee [Hj], suggesting that the UHE cosmic- 
rays with energy from 6 x 10 17 eV above are all probably 
extra-galactic [g, [ll[ . In this paper we show that the 
energetics and spectrum of CRs from the second knee to 
10 19 eV may be due to extra-galactic hypernovae, similar 
to the hyper-energetic supernova SN1998bw. 

SN1998bw was striking not only in its unusually large 
explosion energy, £~ 3 - 5 x 10 52 erg (so called "hyper- 
novae" [IH), but also in that it was associated with a very 



sub-energetic GRB, GRB980425, with an isotropic equiv- 
alent gamma-ray energy E 1 ~ 10 48 erg [l5[. The possible 
transition from extreme SNe to GRBs was implied in the 
magneto-rotational mechanism of the SN explosions pro- 
posed in Ref . (lfjl , \v\ . Now the connection between GRB 
and SNe has been established through observations. The 
hypothesis for SN-associated GRBs such as GRB980425 
that their weakness is an apparent effect due to seeing 
the explosion off-axis runs counter to observational tests 
based on long-time radio observations [l8(, so it is likely 
that this GRB is inherently dimmer than typical. The 
observations of the radio afterglow of this event showed 
that about 10 50 erg of kinetic energy were released in the 
form of a mildly relativistic ejecta Due to the large 
supernova explosion energy and the much lower than typ- 
ical GRB energy, attempts have been made to ascribe 
the GRB event to the shock from the mildly relativistic 
ejecta as it breaks out through the hypernova progeni- 
tor's outer envelope 0, Hlj. A generally accepted con- 
clusion, however, has not yet been reached. A recently 
detected strong thermal X-ray emission component in an- 
other sub-energetic burst (GRB060218), associated with 
SN2006aj , may also be associated with a semi-relativistic 
supernova shock breakout, in which the mildly relativistic 
supernova ejecta has an energy > 10 49 erg [23J. The ra- 
dio observations of this burst, as well as those of another 
hypernova burst, GRB031203/SN20031w, also indicate 
that there is a significant energy in the mildly relativistic 
ejecta j24(. We will use the term semi-relativistic hyper- 
novae to denote such supernovae exhibiting a mildly rela- 
tivistic ejecta component, seen in association with GRBs. 
The recently discovered SN2006gy [25[ also has a large 
explosion energy, but continued multi-wavelength mon- 
itoring has not yielded any evidence for an associated 
GRB, so this object may be in the class of normal hy- 
pernovae with large explosion energy, but without any 
significant mildly-relativistic ejecta [261 ] . 
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FIG. 1: The expected spectrum of CRs as a function of en- 
ergy e produced by hypernova remnants with a distribution of 
the ejecta kinetic energy with velocity Ek oc (F/S) - ". Dashed 
lines indicate the flat, injection spectrum (£ 2 dn/de oc e°) from 
single velocity ejecta for three different velocities. The con- 
volved contribution from different velocity ejecta, denoted by 
the solid line, leads to a spectrum e 2 dn/de oc e~ a ^ 2 . The 
inset shows the kinetic energy distribution of three nearby 
hypernovae associated with sub-energetic GRBs. The data 
points are from Ref . [27l] . The solid and blank circles de- 
note the energy of the slowest ejecta and the mildly rel- 
ativistic ejecta, respectively. A kinetic energy distribution 
Ek oc (T[3)~ a (the solid line) gives a rough fit to the data 
of SN1998bw/GRB9802425 and SN20031w/GRB031203 with 
a ~ 2.4, while for SN2006aj/GRB060218, the slope is a ~ 
-1.7. 



In Figure 1 (inset) we show the kinetic energy dis- 
tribution of the supernova ejecta associated with these 
three nearby sub-energetic GRBs, ranging from the bulk 
of the ejecta at T/3 ~ 0.1 to the mildly relativistic ejecta 
(r/3 ~ 1), where (3 = v/c and T are the ejecta normalized 
velocity and bulk Lorentz factor, respectively. Surpris- 
ingly, even though the energy estimates of the high veloc- 
ity ejecta from the radio observations are crude, all three 
hypernovae give a roughly consistent extrapolation of the 
slope of the kinetic energy distribution from the low to 
the high velocity end. SN1998bw and SN20031w give a 
consistent slope of about ~ —2.4, while SN2006aj gives 
a slightly shallower distribution with a slope of ~ —1.7. 
Note that if the explosion is aspherical, the kinetic en ergy 
released in SN1998bw may be lower, ~ 2 x 10 52 erg [47j . 
which will lead to a slightly shallower slope. It has been 
shown in [23] that the relatively shallow decay of the ra- 
dio afterglow of GRB060218/SN2006aj can be modelled 
with a shock expansion r oc t 0,85 , appropriate for a core- 
collapse supernova explosion with a continuous distribu- 
tion of ejecta velocities 28], propagating into a stellar 
wind environment of density p oc r~ 2 . This provides 
a plausible scenario for a continuous distribution of the 



ejecta kinetic energy over velocities, ranging from the low 
velocity (0.1c) supernova bulk ejecta to the mildly rela- 
tivistic ejecta within the same explosion, of the form 

E k oc (Tl3)- a . (1) 

Such a distribution of velocities is naturally expected in 
an outflow spreading over a wide range of angles [29] . 

Standard hydrodynamic collapse calculations involving 
non-relativistic shocks result in a kinetic energy profile 
Ek oc (r/3) -5 [2l|, with a negligible fraction of the kinetic 
energy at mildly relativistic velocities, consistent with 
the radio observations of, e.g., local type Ib/c supernovae 
19941 and 2002ap. This very steep velocity profile implies 
negligible contribution to the highest energy CRs by high 
velocity ejecta [22j. On the other hand, ultra-relativistic 
shocks result in a much flatter profile, Ek oc (T(3)~ a with 
a ~ 1 [21]. For shocks in the trans-relativistic velocity 
regime, the energy distribution has not been calculated, 
but the above assumed slope of a ~ 2 seems to be inter- 
mediate between the two extreme regimes. An important 
implication of such a continuous energy distribution in 
the semi-relativistic regime is that there is a significant 
amount of energy in the high-velocity ejecta of a hyper- 
nova. At this high velocity, the hypernova blast wave 
could accelerate CRs to energies as high as 10 19 eV, as 
we show below. 

Maximum energy of accelerated CRs. Diffusive 
shock acceleration in supernova shock fronts has been 
extensively studied. The maximum energy for CR ac- 
celeration by the SNRs is usually thought to be limited 
to ~ 10 14 — 10 15 eV for the case of typical interstellar 
magnetic fields of a few ^G. It was suggested in Refs.Q 
that the stellar wind from Wolf-Rayet (WR) star may 
have a relatively high magnetic field with a dominant 
component transverse to the shock normal, so that the 
SN explosion in these winds can accelerate particles to 
a much larger maximum energy. Recently there are also 
suggestions (e.g. [30|) that the magnetic field can be am- 
plified non-linearly through MHD turbulence excited by 
the CRs in the vicinity of the shock to many times the 
pre-shock values, thus significantly increasing the accel- 
eration rate and hence increasing the maximum energy. 
Amplification of the magnetic field to an equipartition 
value is generally assumed in radio SNRs and in GRB 
afterglow shocks, and has gained support from recent 
X-ray observations of several young SNRs[31]. We here 
consider a semi-relativistic hypernova ejecta with a ve- 
locity distribution given by Eq.(l), expanding in the stel- 
lar wind characteristic of WR stars, which are thought 
to be the progenitors of the hypernovae associated with 
GRBs. Different from Refs.[3(, we here consider a ran- 
dom, CR-amplified magnetic field with a strength close 
to the equipartition value. During the free expansion 
phase, the magnetic field is B 2 /8tt — 2tBPw{R)c 2 (i 2 , 
where eg = Q.leg _i is the fraction of the equipartition 
value of the magnetic field energy and p w is the mass 
density of the stellar wind at radius R. The magnetic 



3 



field at the free-expansion radius R is 
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where M is the wind mass loss rate, whose average value 
is 3 x 10~ 5 M Q yr _1 for WR stars, and v w = 10 3 iv, i3 km/s 
is the wind velocity. The maximum energy that can be 
attained in a shock depends on the magnetic field config- 
uration. If the magnetic field is mostly perpendicular to 
the shock surface, the scattering coefficient for particles 
is smaller than in the case of a parallel shock and as a 
result, the maximum energy is given by ~ ZeBR [32| . 
In our case, however, we assume that the CR-amplified 
magnetic field is random in direction and the maximum 
energy of the accelerated particles is (e.g. [3fl l33l|) 
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The stellar wind of WR stars is composed largely of H 
and He (Z = 1 and Z = 2, respectively). In this shock 
acceleration scenario, the maximum energy e max is pro- 
portional to the square of the shock velocity, so a higher 
velocity hypernova ejecta can lead to a higher e max . Note 
also that during the free expansion phase of the ejecta, 
£max is independent of the radius. For the assumed ve- 
locity distribution of Eq.(l), the bulk of the ejecta has a 
velocity of 0.1c and the maximum CR energy correspond- 
ing to this (low-end) velocity ejecta is about 10 17 5 Z eV 
for typical parameters of the stellar wind. 

The spectrum of the CRs. - For a single veloc- 
ity ejecta, the differential spectrum of the accelerated 
protons is given by the injection spectrum, which is 
dN/de oc e 1 with 7 ~ 2.0, for both non-relativistic 
shocks and semi-relativistic shocks [34| . However, if the 
hypernova produces a kinetic energy distribution spread 
over the velocity with the same explosion, as described 
by Eq.(l), the final CR spectrum detected at Earth is 
determined by the kinetic energy distribution profile, as 
illustrated in Fig.l. This can be understood as higher 
energy CRs being contributed dominantly by higher ve- 
locity ejecta, which are represented with a smaller total 
energy. As the maximum CR energy for a particular ve- 
locity is £ max oc (r/3) 2 and the energy distribution of the 

ejecta is Ek cx (r/3)~ Q , we see that Ek oc £ m "i 2 . Con- 
volving the contribution from the different velocity ejecta 
(Fig.l), we expect a final differential energy spectrum of 
the CRs of the form 



''(dN/de) oc e~ a/2 , 



(4) 



where e is the energy of cosmic ray particles, which re- 
lates to the ejecta kinetic energy as E^(e) oc e 2 (dN/de). 
Fits to the observed CR data give a differential spec- 
trum J = C(E/6.3 x loiSeV)- 3 ' 20 "- 05 for 4 x 10 17 eV < 
E < 6.3 x 10 18 eV and J = C(E/6.3 x io 18 e V)- 2 - 75±a2 



for 6.3 x 10 18 eV <£<4x 10 19 eV with C = (9.23 ± 
0.65) x 10 _33 m _2 s _1 sr _1 eV _1 f{|. The observed spectral 
slope at energies 4 x 10 17 eV < E < 6.3 x 10 18 eV implies 
a = 2.40±0.1 in Eq. which is roughly consistent with 
the theoretical slope deduced in Fig. 1 . Note that the en- 
ergy losses due to photo-pair production with CMB may 
become important above ~ 10 18 eV and thus steepening 
the injection spectrum Q. However, the details depend 
on the specific value of the injection spectrum and source 
evolution. In all cases this is a small correction which can 
be accommodated by using a slightly smaller value for a 
and thus a slightly harder injection spectrum from the 
HNe. 

Hypernova rates and the observed flux — We esti- 
mate now how many SN1998bw-like hypernovae per unit 
volume per unit time are needed to produce the CR 
flux from 6 x 10 17 to 10 19 eV. Assuming that the ki- 
netic energy output from one SN1998bw-like hypernova 
is -Efc,HN = 5 x 10 52 erg , the local kinetic energy release 
rate by hypernovae is 



e'k(z — 0) — i?HN-Efc,HN 

= 2.5 x 10 46 ( ^ 



500Gpc-3yr- 



erg Mpc 3 yr 



(5) 



Adopting an efficiency factor 1/6 for the conversion 
of ejecta kinetic energy into CR energy [6J, and 
l/ln(e max /e m i n ) ~ 0.1 as the fraction of the total CR 
energy that is contributed by each decade of energy, 
the local CR energy generation rate per energy decade 
at 10 17 5 Z eV, corresponding to v = 0.1c, is €cr,o — 
0.016e'fc(z = 0). The corresponding expected CR flux is 



e 2 J= (c/AirH )ecRfifz 



(6) 



where Hq = 70 km s 1 Mpc 1 is the Hubble constant 
and 



fz 



Hn 



dz(dt/dz)S(z)(l + z)- 



(7) 



is the correction factor for the contribution from high- 
redshift sources. Here (dt/dz) — H^ 1 (z)/(1 + z), with 
H(z) being the Hubble parameter at cosmological epoch 
z, and z max is the maximum redshift corresponding to 
the mean free path against photopion production in the 
CMB, whose value is z max > 4 for protons with energies 
< 10 19 eV. The value of f z is ~ 2 - 3 for z max = 4 and 
for the source evolution function S(z) given by different 
star formation rates (SFR) in Ref. [38j or the broken 
power-law estimate in Ref. [39| for a standard ACDM 
cosmology. Here we assume f z s» 3. At e = 10 17 - 5 Z eV, 
we get a CR flux of 



j = io- 28 W R ™ 



£ ) eV- 1 m- 2 s- 1 sr- 1 



VSOOGpc-Syr- 1 J \ 3 

(8) 

Comparing this to the observed CR flux of 1.5 x 
10 _28 (e/10 17 ' 5 eV) -3 2 , we infer a required hypernova 
rate of 

i? HN = 750Z- 1 - 2 (/ z /3r 1 Gpc- 3 yr- 1 (9) 
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Assuming Z—l (or 2) and f z = 3, one can derive a 
required hypernova rate of 330 — 750 Gpc _3 yr _1 . Com- 
paring this with the local rate of "normal" type Ib/c SNe, 

- 2 - 5 x 10 4 Gpc^yr- 1 (MUlI, one can find that the 
ratio of the required hypernovae rate to the normal Ib/c 
SNe rate is ~ 1 — 4%, which is consistent with the value 
observed in the local universe ~ 7% [36]. The required 
semi-relativistic hypernova rate is also consistent with 
the observed rate of low- luminosity GRBs [U [H, Hf|. 
Since different SFR give different values of f z , one can in 
principle use the required hypernova rate to constrain f z 
and therefore constrain the SFR. 

Neutrino emission from CRs interacting with hyper- 
nova photons. — The kinetic energy of the highest ve- 
locity ejecta of hypernovae is converted into the highest 
energy CRs in a relatively short time, given by the free 
expansion time t < t^ cc before the ejecta is decelerated 
by the swept-up stellar wind. During this time the hy- 
pernovae remain very bright in the optical-UV band. As 
a consequence, the high-energy protons can interact with 
the hypernova optical-UV photons during the first tens of 
days after the explosion, leading to pry neutrino produc- 
tion. This neutrino production mechanism differs in its 
origins from another, recently suggested neutrino produc- 
tion scheme in low- luminosity GRBs 40] , which involves 
the usual internal shock model for the acceleration of pro- 
tons and production of the target photons. It also differs 
from the suggested mechanism of neutrino production in 
star burst galaxies, where neutrinos are thought to arise 
from pp(pn) interactions between the accelerated protons 
and the interstellar medium [4l[ . 

Around the peak photon luminosity time (typically 
a few tens of days) of type Ib/c supernovae, only the 
high velocity ejecta with T(3 > 1 has been decelerated, 
and only about 2 x 10 50 erg of energy goes into such 
high velocity ejecta. The resulting high-energy protons 
are above the photopion interaction threshold with ther- 
mal hypernova photons if their energy satisfies e p £ 7 > 
0.3GeV 2 , i.e. e p > 10 18 eV for UV-optical photons. 
The cooling rate of a proton of energy e p due to pion 
production via A + resonance is t^} = E^ 1 {dE/dt) ~ 
n(£ 7 )crp eafc c£peak(A£/£ peafe ). Here n(e 7 ) is the number 
density of hypernova photons at the peak of the black- 
body distribution, a pea k = 5 x 10~ 28 cm~ 2 is the cross 
section for pion production for a photon with energy 
e = £peak = 0.3 GeV in the proton rest frame, ^ pea k = 0.2 
is the average fraction of energy lost to the pion, and 
Ae = 0.2GeV is the peak width. Thus, the fraction of 
energy lost by protons to pions is 

fpj — R/(r{3ct P j) {-\r\\ 
= 0.2L SNi 4 3 (i?/10 16 cm)- 1 ^/! eV)- 1 [ ' 

where Lsn — 10 43 erg s _1 is the bolometric luminosity of 
SN1998bw around its peak brightness (see e.g. [l5[) and 
e 7 = 1 eV is taken as the optical photon energy. The v p 
(pn) from pion decays will have a typical energy 0.05e p 



and the expected diffuse neutrino flux is 

- |/P7(4 j p)(r/3>1) = 0-3 x 10- 10 
x ( sooofe-Vi ) (fe) (£) GeV cm^s-isr- 1 

For the thermal spectrum of the target photons, only 
a small energy range of protons with e p > 10 18 eV can 
interact effectively with the photons, so we expect the 
neutrino diffuse emission to peak at ~ 5 x 10 16 eV. The 
maximum neutrino energy would be 1.5 x 10 18 eV corre- 
sponding to the maximum CR energy of 3 x 10 19 eV. The 
probability for a high-energy muon neutrino to interact in 
ice/water producing a high-energy muon detected by the 
embedded instruments is P(ev) ~ 4.6 x 10 _4 (e 1/ /PeV) - 55 
in the PeV-EeV range [42| ■ Integrating the flux in 
Eq. ( 1 1 ) times the probability inthe5xl0 17 -1.5xl0 18 eV 
range, we get an event rate of w 10 -2 yr _1 km -2 over 
2tt sr in a neutrino telescope. Including muon anti- 
neutrinos would increase the rate by a factor 2. This 
rate is too low for cubic kilometer detectors, but future 
> 100 km 2 telescopes such as ANITA and ARIANNA 
may detect these ultra-high energy neutrinos [43| . 

Discussion. — We have proposed that cosmic rays from 
the second knee to 10 19 eV, whose origin has been de- 
bated, may be produced by extra-galactic hypernovae 
similar to SN1998bw, which are associated with under- 
energetic GRBs. The CRs below the second knee may be 
due to heavy ions accelerated in Galactic sources, such as 
Galactic supernovae expanding into the ISM Q or stellar 
winds Q , and Galactic trans-relativistic hypernovae [44| • 
This is supported by the measurements from KASCADE 
[3] that above the "knee" at 3 x 10 15 eV the chemical 
composition is increasingly richer in heavy nuclei. Re- 
cent HiRes measurements [l3[ found that the chemical 
composition changes back towards lighter (proton) com- 
position at and above the "second knee", suggesting a 
transition from Galactic CRs to extra-galactic CRs that 
consist primarily of protons. A smooth transition in the 
CR spectrum between Galactic and extra-galactic com- 
ponents may be reasonable, according to the numerical 
estimates by Hillas (see Fig. 8 of [45|), who found that 
the rapidly falling spectrum may accommodate a factor 
of 3 or more uncertainty in the extra-galactic compo- 
nent without providing a visible clue to the joint point 
of the overall spectrum. Berezinsky et al. @ also consid- 
ered an extra-galactic origin of the cosmic-rays above the 
second knee, similar to us, but with a single component 
extending up to the highest energy that originates from 
AGNs. The so-called "dip" between the second knee and 
the ankle, i.e. first steepening and then hardening of 
the spectrum, has been suggested as being due to adi- 
abatic and pair-production energy losses by cosmic-ray 
protons. In a multi-component model such as ours, the 
"dip" can arise due to a steep injection spectrum, from 
HNe above the second knee (modulo small adiabatic and 
pair-production energy losses) which is taken over by a 
harder spectrum from normal GRBs above the ankle. 
Since hypernovae/under-energetic GRBs belong to the 
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same class as normal GRBs, one should also expect a 
natural flux matching at the ankle from these two com- 
ponents. 

We have shown that semi-relativistic hypernovae can 
accelerate CRs to energies < fO 19 eV, and provide the 
right flux density between the second knee at 6 x 10 17 
eV and 10 19 eV. The assumed ejecta velocity distribution 
profile is consistent with current supernova-GRB obser- 
vations. Confirmation of its theoretical plausibility would 
require further detailed numerical investigations of hy- 
pernova explosions and shock propagation through the 
progenitor envelope. 

Recently, there has been evidence indicating that both 
normal GRBs and sub-energetic GRBs associated with 
hypernovae are preferentially found in low-metallicity 
galaxies [ir| . This would imply that hypernova rates in 
normal metallicity galaxies such as our Milky Way may 
be low, so that the Galactic contribution to the high- 



est energy CRs would be unimportant (48|. Most of 
the flux would thus be expected to originate in distant, 
low-metallicity galaxies and the distribution of these cos- 
mic rays should be isotropic. Since the propagation of 
the cosmic rays in the energy range 6 x 10 17 — 10 19 eV 
may be significantly deflected by the magnetic field in our 
Galaxy, and possibly by the intergalactic magnetic field 
as well, there would be no expected correlation between 
the arrival direction of cosmic rays in this energy range 
and the hypernova host galaxies. 
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